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Abstract 

We study the effects of decoherence on the entanglement generated by 
Unruh effect in accelerated frames by using various combinations of an 
amplitude damping channel, a phase damping channel and a depolariz- 
ing channel in the form of multilocal and collective environments. Using 
concurrence as entanglement quantifier, we show that the occurrence of 
entanglement sudden death (ESD) depends on different combinations of 
the channels. The ESD can be avoided under a particular configuration 
of the channels. We show that the channels can be used to distinguish 
between a moving and a stationary frame. 
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1 Introduction 

The last two decades has witnessed the rapid development in the field of quan- 
tum information and strengthened the notion that entanglement is not only 
one of the fundamental concepts of quantum theory but also a key concept for 
transmitting and processing quantum information [T]. Entanglement between 
spatially separated parties is used as a potential source for quantum teleporta- 
tion of unknown states ^ , quantum key distribution 3 , quantum cryptography 
[4] and quantum computation [5l[6]. However, the behavior of entanglement be- 
tween various systems is still not fully known and efforts are on the line to un- 
derstand its dynamics deeper under various setups. The study of entanglement 
in a bipartite system has recently been extended to the relativistic setup and 
its behavior for various fields has been examined [71 |H1 [HI UHl HU [H]- However, 
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these investigations on entanglement are limited to closed quantum systems. 
Practically, quantum systems are influenced by its environment that may re- 
sults in non-unitary dynamics of the system. For most accurate and practical 
results in quantum processing, the effect of environment on the entanglement 
between spatially separated systems needs to be necessarily investigated. The 
environmental effect on a quantum system gives rise to the phenomenon of de- 
coherence that causes an irreversible transfer of information from the system to 
the environment [T31 [TS] . 

Alsing et al [7] have shown that the entanglement between two modes of 
a free Dirac field is degraded by the Unruh effect and asymptotically reaches 
a nonvanishing minimum value in the infinite acceleration. In Ref. jl6) , the 
dynamics of entanglement for Dirac field in the presence of Unruh effect has 
been studied under amplitude damping channel. The effect of decoherence on 
the entanglement of Dirac field using a phase damping channel, a depolarizing 
channel and a phase flip channel has been investigated in Ref. [TT] . 

In this paper we investigate the effect of decoherence on the entanglement 
of Dirac field in a noninertial system by considering various combinations of an 
amplitude damping, a phase damping and a depolarizing channel in the form 
of multilocal, collective and global noises. For example, we consider that if one 
qubit is locally coupled to an amplitude damping channel and the other is locally 
coupled to a phase damping channel then both the qubits are collectively coupled 
to a depolarizing channel. Using concurrence as the entanglement quantifier, 
we show that the rate of decrease of concurrence depends on the coupling of 
the kind of a channel to a particular qubit. For example, under the influence of 
multilocal coupling when Alice's qubit is influenced by phase damping channel 
and Rob's qubit is influenced by amplitude damping channel, the entanglement 
sudden death (ESD) [TS] occurs earlier for large acceleration of the frames under 
some particular conditions. On the other hand. Under some other setups, the 
ESD can be avoided or delayed or even faster. 

We rigorously consider three different cases. In one case we allow Rob's 
qubit to interact locally with an amplitude damping channel and Alice's qubit 
interact locally with a phase damping channel. In the second case, the coupling 
of these two channels are reversed in the sense that Rob's qubit interacts with a 
phase damping channel and Alice's qubit interacts with an amplitude damping 
channel. However, in both these cases the two qubits are coupled collectively 
with a depolarizing channel. In the third case, the two channels in the multi- 
local coupling of the qubits are depolarizing and phase damping channels while 
collectively the two qubits are coupled to an amplitude damping channel. 

We consider that Alice and Rob, the two parties, share the following maxi- 
mally entangled state at a point in flat Minkowski spacetime 

= ^ {\Ou^a)m\O^Jm + \Iu:a)m\Iu:Jm) ■ (1) 

In Eq. ([T]) |0) and |1) kets represent, respectively, the vacuum and excited states. 
The first entry in each pair of the kets is in Alice's possession and the second 
entry is in Rob's possession. The subscirpts lun {N ~ A, R) of the kets reflect 
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that we consider the initial entanglement between these two modes and all the 
rest modes are in vacuum states. Furthermore, we suppose that each player has 
a device that is capabale of decting only the mode that she/he has shared for the 
initial entanglement. Rob then moves with a constant acceleration and Alice 
stays stationary. For accelerated observer the suitable coordinates are Rindler 
coordinates. These coordinates define two different regions, usually denoted by 
/ and //, and are called Rindler regions. An observer in one region is causally 
disconnected from the other region. In other words, an observer in region / 
has no access to information that leaks out to region // and vice versa (for 
detail see [7] and reference therein). A given Minskowski mode of a particular 
frequency spreads over all positive Rindler frequencies that peaks about the 
Minskowski frequency |19[ 120] . However, to simplify our problem we consider a 
single mode only in the Rindler region /, which is valid if the observers' detectors 
are highly monochromatic that detects the frequency uja ~ = From this 
point onward, with this approximation, the frequency subscript of kets will be 
dropped. 

From the accelerated Rob's frame, the Minkowski vacuum state is found to 
be a two-mode squeezed state [7] , 

|0>Af =cosr|0)/|0)// + sinr|l)/|l)//, (2) 

where cosr — i^^-'^'^^'^l ^ i) . The constant cj, c and a, in the exponential 
stand, respectively, for Dirac particle's frequency, the speed of light in vacuum 
and Rob's acceleration. In Eq. ([2]) the subscripts / and // of the kets repre- 
sent the Rindler modes in region / and //, respectively. The excited state in 
Minkowski spacetime is related to Rindler modes as follow [7] 

|1)a./ = |1>/|0)//. (3) 

In terms of Minkowski modes for Alice and Rindler modes for Rob, the 
maximally entangled initial state of Eq. ([1]) by using Eqs. ^ and ([3]) becomes 

\'^)aj.u = ^ (cosr|0)A|0)/|0)// + sinr|OU|l)/|l>// + |1>a|1>/|0)//) . (4) 

Traditionally, we consider Rob to be in region /, then, he is causally disconnected 
from region //. All the modes in region // of Eq. (U) needs to be discard. So, 
tracing over all the modes in region // leaves the following mixed density matrix 
between Alice and Rob 

PAj = i[cos2r|00)A/(00|-f cosr(|00)A,/(ll| + |ll)A/(00|) 

sin2r|01)^j(01| + |ll)^,/(ll|]. (5) 

2 The system in noisy environment 

The effect of decoherence is studied via quantum channels in the Kraus operators 
formalism. The Kraus operators for the three channels of a single qubit system 
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Table 1: A single qubit Kraus operators for amplitude damping channel, phase 
damping channel and depolarizing channel. 
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that we use in this paper are given in Table 1. When a single qubit system 
evolves under the action of amplitude damping channel there is a probability 
to change the state of the system from state |1) to state |0). When the density 
matrix of a single qubit system is influenced by a phase damping channel, the 
diagonal elements remains unaffected while the off-diagonal elements decay. On 
the other hand, the study of the effect of a depolarizing noise on a quantum 
system is important because a pure state of a system when influenced by such 
a noise is replaced, with some probability, with a maximally mixed state. In 
Table 1, the action of the Kraus operators of each channel is parameterized by 
the decoherence parameters pi {i — 1,2,3) whose values lie between and 1. 
For Pi = 0, the channels have no effect on the system and for = 1, the system 
is fully decohered. Furthermore, the Kraus operators for each channel obey the 
completeness relation ElEi = I. The evolution of the initial density matrix 
of the system when it is influenced by the global environment is given as follow 

/'Z - E E E {Et^EfE^) p,, {Et'EfEt'''^) , (6) 

i j k 

where E^ = E^^h, Ef = h^E^ are the Kraus operators of the local coupling 
of Alice's qubit and Rob's qubit, respectively. The I2 is identity operator and 
the subscripts m and n represent the number of single qubit Kraus operators 
of the particular channel that is coupled to Alice's qubit or Rob's qubit. The 
Ef^ are the Kraus operators for the case of collective coupling and are formed 
from all the possible combinations of the tensor product of the Kraus operators 
of a single qubit channel in the form E^ (g) E^. The subscripts g, represents the 
number of Kraus operators for a one qubit channel used as a collective environ- 
ment. The system is said to be coupled with multifocal environment when both 
the qubits are influenced by their local environments and it is said to be under 
the action of global environment when both multifocal and collective couplings 
are switched on simultaneously (global — multifocal + collective). First we con- 
sider the situation in which Alice qubit is locally influenced by phase damping 
channel, Rob's qubit is locally influenced by amplitude damping channel and 
the two qubits are then coupled collectively to depolarizing channel. Under such 
condition, the summations over i, j and k in Eq. Q are, respectively, given by 
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i — 0,1, 2, ...16 and j, fc = 0, 1. If pi, p2 and stand for decoherence param- 
eters of amplitude damping channel, phase damping channel and depolarizing 
channel, respectively, then the non-zero elements of the final density matrix of 
the system in the light of Eq. ^ become 



Pii = ^{-SpUpi - 1) - (2P3(4P3 -9)+ 9)(pi - 1) cos2r 
-6p3{pi + 1) + 9(pi + 1)), 

Pi4 = ^(3 - 4p3) V(-l +P2) cosr. 



P22 = ^{9 + 8pl{-l+pi)~9pi+6p3{l+pi) 



1 

36 ~ 

+ (9 + 2p3{-9 + 4p3))(-l cos2r), 
P33 = 7^(P3(9 + 4p3(-l - 15pi) + 9pi 

+P3{-3 + 4p3)(-l + Pi) COS 2r), 
P41 - ^(3-4p3)V(-l+Pi)(-l+P2)cos2r, 

P44 = 4(9-9pi+p3(-9-4p3(-l+pi)-M5pi) 

io 

-p3(-3-}-4p3)(-l+Pi)cos2r). (7) 

The spin- flip matrix of the final density matrix of Eq. ([6]) is defined as pf = 
((72 ® (^2) Pf {cf2 ® 02), where 02 is the Pauli matrix. The degree of entangle- 
ment in the two qubits mixed state in a noisy environment can be quantified 
conveniently by concurrence C, which is given as [211 122) 

C = max|o, /aT- v/A^- /Xl- v/aI} A, > A,+i > 0, (8) 

where are the eigenvalues of the matrix p//5/. The eigenvalues under the 
action of global environment become 

Ai,2 = ^ (54p3 - i&pl + 81pi - 216p3Pi + 144p^pi - Sip? + 216p3P? 

-IUpIpI + (3 - 4p3)'(-l +Pi)(-24p3(-l +P2) + lQpl{-l+P2) 
+9(-2 pi + P2)) cof? r + 2(3 - 4p3)2p3(~3 + 2p3)(-l + pif cos^ r 
±2[(3 - 4p3)*(-l +Pi)(-1 + P2) cos^ r(-9(-6p3(l - "^Pif 
+4p2(i _ 2p^f + 9(-l +pi)pi) + 9(3 - Ap^fi-l+Pif cos^ r 
+2(3 - 4p3)'p3(-3 + 2p3)(-l cos'r))]"\ 
A3 = A4-^(3(-96pi(-l+pi)2-H32pt(-l+pi)2-54(-l+pi)pi 

-6p2(_7+ 14pi +pI) + 9p3(5 - lOpi + 13p2)) + 2(3 - Ap^f 

(-l+pi)(-6p3(-l+pi) +4p§(-l+pi) +9pi)cos2r 

+(3 - 4p3)2p3(-3 + 2p3)(-l + pif cos4r), (9) 
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where in Ai^2, the two eigenvalues are differentiated by ± sign, respectively. Us- 
ing Eq. ([8]) the concurrence can straightforwardly be calculated. It is important 
to note that the concurrence defined by the eigenvalues of Eq. © reduces to 
the result of Ref. ^ when the decoherence parameters pi = 0. 




(a) 



(b) 

Figure 1: (color online) (a) The concurrence for different values of the acceler- 
ation is plotted against the decoherence parameter p2, for the case when both 
the qubits are coupled to their respective environments in multifocal coupling. 
The other parameters are set to pi ~ 0.5, and p^ ~ 0. (b) The concurrence 
against pi for p2 — 0.5 and ps — for different values of the acceleration of 
Rob's frame. 

To see how the concurrence and hence the entanglement change when the 
system is coupled to noisy environments in the presence of Unruh effect, we plot 
the concurrence for each kind of the coupling of the three channels against the 
decoherence parameters for various values of the acceleration r of the noninertial 
frame. In Fig. l(a, 6), the concurrence, when the collective environment is 
switched off (ps = 0), is plotted for the case of multifocal coupling. Fig. la 
shows the behavior of concurrence against the decoherence parameter p2 (phase 
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damping) for pi = 0.5 (amplitude damping). One can see that the concurrence is 
strongly acceleration dependent as the initial value of the concurrence decreases 
with increased acceleration. For a given acceleration, the decoherence parameter 
considerably damps the concurrence as the system evolves and it becomes zero 
at a definite value of the decoherence parameter (particular time) and hence 
loss of the entanglement occurs. The BSD occurs earlier for the case of larger 
acceleration. On the other hand, the behavior of concurrence is different when 
it is plotted against pi for a particular value of p2- Fig. 16 shows such a 
situation for p2 = 0.5. Here the concurrence is not as much strongly acceleration 
dependent as in the case of Fig. la, except for the larger values of acceleration. 
For smaller acceleration of the noninertial frame, the BSD does not happen, 
the concurrence goes to zero only when the channel is fully decohered (infinite 
time). However, for large values of the acceleration (r — 7r/4), the concurrence 
is influenced identically by both pi and p2- From the different behavior of 
concurrence and hence of BSD in the range of smaller acceleration, the two 
channels can be used to identify the stationary and the accelerated observers. 
The behavior of concurrence is shown in Fig. 2 (a, b) for the case when both 
the collective and multilocal environments are simultaneously switched on, that 
is, when the system is coupled to global environment. In Fig. 2a, the values of 
decoherence parameters that parameterize the multilocal environment are set 
to pi = p2 = 0.1. It can be seen that the concurrence is strongly dependent 
both on acceleration and decoherence parameter. The BSD occurs faster than 
in the case when the system is coupled to multilocal environment only. Fig. 2b 
shows the situation in which pi = P2 = P3 = P- One can see from the figure that 
under such constraint, the concurrence, as compared to Fig. 2a, becomes more 
dependent on decoherence parameter and less dependent on the acceleration. 

Now, we want to see the effect of interchanging multilocal channels that is, 
when Alice's qubit is locally coupled to an amplitude damping environment and 
Rob's qubit interacts locally with a phase damping environment. For this case, 
the summations over i,j and k in Bq. ^ remain unchanged. The concurrence 
can straightforwardly be calculated as in the previous case. Without writing 
its mathematical form, we prefer to show its behavior by plotting it against the 
decoherence parameters. The concurrence as a function of pi for p2 — 0.5 is 
shown in Fig. 3 for different values of acceleration when the collective envi- 
ronment is switched off. The figure shows that the concurrence is significantly 
both acceleration and decoherence parameter dependent such that it decreases 
with the increasing values of both of these quantities. However, no BSD occurs 
for any value of the acceleration, the concurrence goes to zero for every accel- 
eration only when pi — 1. Unlike the previous case, the plot of concurrence 
against P2 is exactly the same as given in Fig. 3. This means that the rapid 
fall of concurrence in the previous case (Fig. la) is somehow compensated by 
the acceleration of the moving frame in this case. The two channels are not 
differentiable in this case. Nevertheless, BSD does occur in this case when the 
system evolves in the global environment. The effect of global environment on 
the concurrence is shown in Fig. 4. A comparison with Fig. 2a shows that the 
concurrence is influenced almost identically in the two cases by the decoherence 
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parameters, however, it is less affected by the acceleration in the later case. 

Next, we consider the setup in which Alice's qubit interacts locally with a 
depolarizing channel and Rob's qubit interacts locally with a phase damping 
channel as before. The two qubits then interact collectively with an amplitude 
damping channel. The summations over i, j and k in Eq. (jS]) become i = k = 
0,1,2,3 and j = 0,1. Using Eq. © the final density matrix can easily be 
calculated. After finding the spin-flip matrix, the four eigenvalues of the matrix 
pp are given by 

Ai.2 = 4(-l+]5i)2(18 + 9p?-9p2-36p3 + 12piP3 
So 

-Upjps + 24p2P3 + 20pl - Spipl + 4pIpI - I6P2P3 
±V2[-(-l+j32)(3-4p3)2cos2r(-3+j33 +P3 cos2r) 
X (-3(1 +PI + 2p\) + 2(-l + pifp^ + (-1 + pi) 
X (3 + 2(-l + pi)p3) cos 2r)]i/2 + (9(_2 +pi+ P2) 
+6(7 - 3pi + 2pl - 4p2)P3 - 8(3 - 2pi + _ 2p2)pl) sin^ r 
+2(-l + pi)p3(3 + 2(-l +pi)p3) sin^ r), 

A3 = A4 = ^(-l+pi)2(3(2p2(_i+pi)2 + 6pi(l + 2pi) 
+p{l + 7pi - 8pl)) + 2(4p2(_i +pi)2 _ gp^ 

-12p(-l + cos2r + p3(3 + 2p3(-l + Pi))(-1 + Pi) cos 470.0) 

The concurrence can be found by using these eigenvalues in Eq. To 
study the influence of the noisy environment on the concurrence, we proceed by 
plotting it against the decoherence parameters as done in the previous cases. 
First we consider the effects of decoherence in the multifocal coupling. The plot 
of concurrence against ps {p2 = 0.2) and p2 {ps = 0.2) for pi = is shown in 
Fig. 5(a,6). Though the behavior of concurrence is different in both of these 
figures, however, the ESD cannot be avoided. The ESD in this case occurs 
faster as compared to the previous two cases. On the other hand, the damping 
of concurrence due to acceleration is stronger in Fig. 56 than in Fig. 5a. Finally 
we show the behavior of concurrence when the system is coupled with the global 
environment. Fig. 6 shows the concurrence of the system against pi for different 
values of acceleration with P2 = P3 = 0.2. It is obvious from the last two figures 
that when one qubit is locally coupled to depolarizing channel there is no way 
to avoid ESD and it happens faster than the other two cases. 

3 Conclusion 

In conclusion, we have studied the behavior of entanglemnent of Dirac fields 
by considering different combinations of the three different quantum channels. 
We show that under multifocal environment when Alice's qubit is coupled to a 
phase damping channel and Rob's qubit is coupled to an amplitude damping 
channel, the ESD may or may not occur. However, the distinguishing effects 
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of the two channels may be used to identify which of the observer is in mo- 
tion. On the other hand, if the multilocal couphng of the two channels with two 
qubits is interchanged, the two channels affect the concurrence identically and 
cannot be distinguished from each other. However, the ESD can be avoided. 
The concurrence goes to zero only when one of the channels is fully decohered. 
Furthermore, it is shown that irrespective of the multilocal coupling's configu- 
ration, the ESD occurs when the collective environment is present. Moreover, 
It is shown that ESD occurs faster and can not be avoided when one of the two 
qubits is locally coupled to the depolarizing channel. 
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(a) 




(b) 

Figure 2: (color online) (a) The concurrence for different values of the acceler- 
ation is plotted against the decoherence parameter p3, for the case when both 
the multifocal and collective environments are switched on. The other param- 
eters are set to pi — P2 — 0.1. (b) The concurrence is plotted for different 
values of the acceleration against decoherence parameter p under the conditions 

Pl=P2=Pd.= P- 
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Figure 3: (color online) The concurrence for different values of the acceleration is 
plotted against the decoherence parameter pi , for the case of multilocal coupling 
when Alice's qubit is coupled to amplitude damping channel and Rob's qubit to 
phase damping channel. The other parameters are set to p2 = 0.5, and ps — 0. 
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Figure 4: (color online) The concurrence for different values of the acceleration is 
plotted against the decoherence parameter , for the case of global environment . 
The other parameters are set to pi = P2 = 0.1. 
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(a) 




(b) 

Figure 5: (color online) (a) The concurrence lor different values ol the acceler- 
ation is plotted against the decoherence parameter p^, for the case when both 
the qubits are coupled to their respective environments in multifocal coupling. 
The other parameters are set to p2 — 0.2, and pi ~ 0. (b) The concurrence 
against p2 for P3 = 0.2 and pi = for different values of the acceleration of 
Rob's frame. 
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Figure 6: (color online) The concurrence for different values of the acceleration is 
plotted against the decoherence parameter pi , for the case of global environment. 
The other parameters are set to P2 = Ps = 0.2. 
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